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Abstract 

We derive a generalized Balitskii-Fadin-Kuraev-Lipatov equation, which apphes directly 
to the perturbative QCD component of total cross section. With the gluon correlation radius 
Rc ~ 0.4fm we reproduce the empirical rate of growth of the hadron-nucleon total cross 
sections. The simultaneous estimate of the triple-pomeron coupling also agrees with the ex- 
periment. 
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Lipatov and his collaborators have shown [1-3] that for short-distance interactions the 
QCD pomeron has the intercept 

aip = 1 + A]p = 1 + ^^^^^^5 , (1) 

TT 

which gives a very large Ajp ~ 1 even with the reasonably small strong coupling as = 
Qs/^T^ ~ 0.4 appropriate for the already short distances (r ~ O.lSfm) (here gs is the color 
charge). On the other hand, the atot oc s^ip^''^) fit of the hadronic total cross sections yields 
A-p{hN) ~ 0.1 [4] (here s is the square of the c.m.s. energy). 

In this paper we derive a particularly simple generalization of the Balitskii-Fadin-Kuraev- 
Lipatov (BFKL) equation directly for the total cross sections. We use the s-channel approach 
to the pomeron, which is based on the technique of light-cone multiparton wave functions 
introduced by the two of the present authors [5,6]. It allows us to easily introduce in the 
gauge-invariant manner the effect of a finite radius for gluon correlations Rc and to evaluate 
the effective intercept of the pomeron for the hadronic scattering regime. We relate the growth 
of the total cross section to the rise of the number of perturbative gluons in the lightcone 
hadrons. The lightcone wave functions of the multiparton states and their interaction cross 
sections were derived in our previous paper [6] and applied to an analysis of the diffractive 
deep inelastic scattering in the Double-Leading-Logarithm Approximation (DLLA) (see also 
[7]). In this paper we extend the considerations of [6] to the BFKL regime, discuss the 
connection between the BFKL and DLLA regimes, and estimate Ajp for hadronic scattering 
processes. 

Our starting point is the lowest order perturbative QCD cross section for the scattering 
of the two color dipoles r and R (here r , R are 2-dimensional vectors in the impact parameter 
plane) 

ao(f, R) = ^ J ^p^^2 y ^s [l - exp(-iA^f)] [l - exp{ikR)^ . (2) 
Here the effective mass of the gluon hg serves as a reminder that the colour forces can not 
propagate beyond the the gluon correlation radius Rc = I /fie and a| must be understood 
as a5(max{/c^, ^})Q;5(max{fc^, -^}) [5]. At r <^ R < Rc, Eq. (0) gives the driving term of 
the DLLA cross section (for a detailed discussion of the DLLA regime see [6]) 

ao{r,R)^Cr^as{r)L{R,r) , (3) 



which is independent of Re- Here L{R,r) ^ log[as{R)/<ys{r)]- In terms of the dipole-dipole 
cross section (0) the perturbative part of the total cross section for the interaction of mesons 
A and B equals 



a^^'\AB) = {{airA,rB))A) 
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j dzAd'^rAdzBd'^rB\'^{zA,rA)\^\'^{zB,rB)\^(To{rA,rB) ■ (4) 

The advantage of the representation (^) is that it makes full use of the exact diagonalization of 
the scattering matrix in the dipole-size representation. Hereafter we discuss a{f, R) averaged 
over the relative orientation of dipoles. 

The perturbative qqg Fock state generated radiatively from the parent colour- singlet qq 
state has the interaction cross section crsi^r, pi, P2) = |[co(pi) + cro(P2)] — |c"o('^)) where pi 2 
are separations of the gluon from the quark and antiquark respectively, P2 = Pi + ^ [6]- The 
increase of the cross section for the presence of gluons equals (we suppress the target variable 
R) 

9 

Aag{r, pi, P2) = (T3{r, pi, P2) - (To{r) = -[ao{pi) + ao{p2) - Mr)] , (5) 

o 

The hghtcone density of soft, Zg 1, gluons in the qqg state derived in [6] equals 
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Here gs{r) is the running colour charge, 2 = niin{r, pi^2}, Ki{x) is the modified Bessel 
function, Zg is a fraction of the (lightcone) momentum of qq pair carried by the gluon, and 
/ dzg/zg = log{s/so) = ^. With allowance for the qqg Fock state the dipole cross section 
takes the form <Jtot{C,, r) = cro(r) + <Ji{r)^, where [6] 

(^i{r) = J d^pi Zg|$i(f,pi,p2,2;g)pAcrg(r,pi,p2) = K ® ao{r) . (7) 

To higher orders in ^, (t(^, r) = J2n=o '^(^ni^)^", where <7n+i = /C ® (t„, so that 

^ = (8) 

is our generalized BFKL equation for the dipole cross section. We emphasize that the in- 
troduction of the gluon correlation length Rc in the kernel /C does not conflict the gauge- 
invariance constraints a{r) and |$(ri, pi, p2, ^ at r — > 0, and A(7g{r, pi, P2) 



at Pi, 2 —>■ 0. The essential ingredient of this derivation is the subtraction of cro(r) in Eq. 
which in a simple and intuitively appealing form takes care of the virtual radiative correc- 
tions. Once the dipole-dipole scattering problem is solved, Eq. (^) gives the hadron-hadron 
cross section. 

In the BFKL scaling limit of r, pi, P2 <^ Rc and fixed as 
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the kernel JC becomes independent of Rc and with the fixed 0:5 it takes on the scale-invariant 
form. The BFKL eigenfunctions of Eq. (H) are E(to,^,r) = (r2)5+'^ exp[^A(ci;)] with the 
eigenvalue (intercept) [here r = rn, pi = rx and p2 = r{x + n)] 



J aP{x + ny n Jo z{l — z) 

= ^-^[2^{i)-nl-^)-nl+^)], (10) 
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where \E'(x) is the digamma-f unction. The final result for A(c<j) coincides with eigenvalues of 
the BFKL equation. Indeed, in the scaling limit of /ic — ^ our Eq. (§) can be transformed 
to the same form as the original BFKL equation for the differential distribution of gluons 
[1,2]. 

When uj is real and varies from — | to and to 1, also the intercept A(co') is real and 
varies from +00 down to A(0) = Aip and back to +00, along the cut from j = 1 + Ajp to 
+00 in the complex angular momentum j plane. \i lo = iv and v varies from —00 to and 
to +CXD, then the intercept A{iu) is again real and varies from — cxd up to A(0) = Ajp and 
back to —00, along the cut from j = — 00 to j = 1 + Ajp in the complex j-plane. The choice 
of the latter cut is appropriate for the Regge asymptotics at ^ ^ 1 and the counterpart of 
the conventional Mellin representation is 

/+00 

/ + 00 
dvf{v) exp[2zz/ log(r)] exp(A(iz/),^) (11) 
-00 

where the spectral amplitude /(z/) is determined by the boundary condition = 0,r): 

f{iy) = i| rfr^^^exp[-2^z/log(r)] . (12) 



In the BFKL regime, the rightmost j-plane singularity corresponds to the asymptotic cross 
section 

o-ip(^,r) oc rexp(^Aip) . (13) 
The solution of Eq. (^ can be written as 

/dv' 
P^K(e,r,r')a(e = 0,r'), (14) 

where in the BFKL regime the evolution kernel equals 



1 /■ r r 

K(^,r,r') = — duexp 2iiy\og — 
Tc J I r' 
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exp(AipO / (logr-logrpn 
" ' n eA"(0) ) ■ ^''^ 

The 'diffusion' kernel (plSf ) makes it obvious that starting with = 0,r) which was con- 
centrated at the perturbative small r < R <^ Rc one ends up at large ^ with a{^, r) which 



extends up to the nonperturbative r ~ i?exp(y ^A"(0)) > Rc . This 'diffusion' towards 
large r is further accelerated if the running coupling is introduced. Thus, the scattering of 
even very small dipoles of size r ~ i? <^ _Rc and/or deep inelastic scattering at ^ Rc^, 
will eventually be dominated by interactions of the perturbative gluons sticking out of small 
dipoles at a distance p ^ Rc, and Aip will be determined by interactions at the scale Rc and 
by the frozen coupling ag'^\ 

The analytic solution of Eq. (8) at finite Rc and with the running coupling is not available, 
and we resort to the numerical analysis. We use the running QCD coupling as{r) = 67r/[(33 — 
2Nf) log(l/Ar)] with A = 0.2GeV, and at large r we impose the simplest freezing as{r) = 



a 



(fr) 



min{a5'(_Rc), 1} • Firstly, we calculate Acff{r, R) = cri(r, R)/(jQ{r, R) for the boundary 



condition (t{C, = 0, r) given by the dipole-dipole cross section @. The first interesting case is 
the scattering of the equal-size dipoles (Fig. 1). Here the kernel /C suggests the quasiclassical 
estimate for the effective intercept 

A^{r) = ^^^as{r). (16) 

TT 

Indeed, we find that at small r the ratio P = Aejj(r, r)/Aip(r) tends to a constant j3 ~ 0.57 
independent of the gluon correlation radius Rc- The effective intercept Ae//(r, r) rises with 
r up to r ~ Rc, then decreases and flattens at r ^ Rc, where both cro,i{r) oc R^. 



Another interesting regime is the DLLA of unequal dipoles r <^ R. In this hmit Eq. (|^ 
takes the form [6] 

a„+i(r) = 1C® fT„(r) = 1^ / -^an{p) , (17) 

which is equivalent to the GLDAP evolution equation [8] and to the first order in ^ gives 
^DLLA{r, R) = ^\og[as{R)/(ys{^)]- Like any logarithmic estimate, this formula works up to 
a constant term ~ 1. In Fig. 2 we show our results for Ae//(r, R) for R = If. The difference 
6 = Ae//(r, R) — ADLLA^r, R) flattens at r/R < 0.2, which is a signal of the onset of DLLA. 

The position of the rightmost singularity in the j-plane can easily be estimated from 
the asymptotic behavior of the numerical solution of Eq. (|^. At R^ = 0.4, 0.28, 0.22 f, 
i.e., at the frozen coupling a^J^^ = 1.0, 82, 0.63, we find estimates Ajp ^ 0.52, 0.41, 0.36, 
respectively. The detailed discussion of the convergence to, and properties of, the limiting 
cross section <7ip{C,,r) will be presented elsewhere. We only notice, that these estimates are 
significantly below the CoUins-Kwiecinski lower bound Ajp > S.Qog'^'^ /n [9]. (The derivation 
of this bound in [9] is flawed by the infrared cutoff which breaks the initial sjnnmetry of the 
BFKL kernel in the momentum space.) 

The two cases of certain theoretical, although of little practical, interest are worth of 
mention. The first is the case of finite Rc at fixed 0:5, the second case is of pc — > in the 
wave function (^ while keeping finite Rc in the strong coupling. Both share the property of 
restoration of the scaling invariance of the kernel /C on the infinite semiaxis logr < log-Rc 
or logr > log-Rc; where as{r) freezes, respectively. On the corresponding semiaxis, the 
eigenfunctions are essentially identical to the BFKL set, the spectrum of eigenvalues will 
evidently be continuous and the j-plane partial waves will have the cut in the j-plane. Here 
we differ from Lipatov [3], who concluded that the running coupling leads to the discret 
spectrum of eigenvalues and to a sequence of poles in the j-plane. The numerical analysis 
shows that the tip of the cut in the j-plane is very close to Ajp as given by Eq. ([I|) with 
as = Os'^^'i ^Le more detailed analysis is needed to check a possibility of a finite departure 
from Eq. (|1|) which depends on the value of al/'^^ 

Finally, let us consider the nN interaction as the typical hadronic scattering. The plausi- 
ble assumption is that the growth of the hadronic cross sections is dominated by the perturba- 
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tive gluons. In this case, af '{nN) = ((ai(r^, rAr))Ar)7r and Aip(7rA^) fa ai {nN) / atotij^N). 
In Fig. 3 we present our prediction for the perturbative QCD contribution to the total cross 
section ao{7rN) and the Ae//(7rA^) vs. the gluon correlation radius Re- We reproduce the 
empirical value Ajp{hN) ~ 0.1 at Rc ~ 0.4fm, when ~ 40% of atot{7^N) is of the perturbative 
origin. 

Besides the rise of the total cross section, variations of the dipole cross section for the 
presence of gluons also contribute to the triple-pomeron coupling A^jp. The method for 
calculating A^jp was presented in [6,7]. With a correlation length Rc ~ 0.4fm, we find 
A^jp{nN) ~ 0.04GeV~^, which is conistent with the experimental determinations [10]. 

In conclusion, we have derived generalized BFKL equation for the total cross sections 
with allowance for the finite gluon correlation radius Rc- We presented the first estimates of 
the perturbative QCD contribution to the rate of the growth of CtotipN) and to the triple- 
pomeron coupling, which are consistent with experiment if Rc ~ 0.4 fm. (Incidentally, the 
instanton model of the QCD vacuum and the lattice QCD calculations give very close value 
of Rc [11].) The irrefutable advantage of having the equation for the total cross section 
and of using the dipole-size representation which diagonalizes the scattering matrix, is that 
they allow an easy incorporation of the unitarity constraints. To this end, we recall that 
consideration of the unitarity corrections in the DLLA limit has already lead to an important 
conclusion [6] that the unitarity correction to structure functions in the diffractive deep 
inelastic scattering at small x satisfies the linear GLDAP evolution equations. 
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Fig.l - The effective intercept Ae//(r, r) for the scattering of two identical dipoles of size r in 
comparison with Aip(r) Eq. (16). The lower box shows the ratio Ae//(r, r)/Aip(r). 
The curves a), b) and c) are for //q = 0.3, 0.5, 0.7 GeV, respectively. 

Fig.2 - The effective intercept Ae//(r, i?) for the scattering of unequal dipoles of size r and 
R in comparison with DLLA formula A.dlla{t^, R)- The curves a), b) are for /iq — 
0.3, O.SGeV, respectively. 

Fig. 3 - The effective intercept Ajp{nN), the perturbative QCD contribution to the total cross 
section ao{TTN) and the effective triple-pomeron coupling for the pion-nucleon scattering 
vs. IIG- 
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